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Abstract 

A precision measurement of the double-differential production cross-section, d^a^^ /dpdfl, for 
pions of positive charge, performed in the HARP experiment is presented. The incident particles are 
protons of 12.9 GeV/c momentum impinging on an aluminium target of 5% nuclear interaction length. 
The measurement of this cross-section has a direct application to the calculation of the neutrino flux of 
the K2K experiment. After cuts, 210 000 secondary tracks reconstructed in the forward spectrometer 
were used in this analysis. The results are given for secondaries within a momentum range from 
0.75 GeV/c to 6.5 GeV/c, and within an angular range from 30 mrad to 210 mrad. The absolute 
normalization was performed using prescaled beam triggers counting protons on target. The overall 
scale of the cross-section is known to better than 6%, while the average point-to-point error is 8.2%. 
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1 Introduction 



The objective of the HARP experiment is a systematic study of hadron production for beam momenta 
from 1.5 GeV/c to 15 GeV/c for a large range of target nuclei 0. The main motivations are: a) to 
measure pion yields for a quantitative design of the proton driver of a future neutrino factory, b) to 
improve substantially the calculation of the atmospheric neutrino flux and c) to provide input for the 
flux calculation of accelerator neutrino experiments, such as K2K and MiniBooNE. 

The measurement described in this paper is of particular relevance in the context of the recent results 
presented by the K2K experiment [21 Ej j which have shown evidence for neutrino oscillations at a con- 
fldence level of four standard deviations. The K2K experiment uses an accelerator-produced beam 
with an average energy of 1.3 GeV directed at the Super-Kamiokande detector. The K2K analysis com- 
pares the observed spectrum in Super-Kamiokande, located at a distance of about 250 km from the 
neutrino source, with the predicted spectrum in the absence of oscillations. This, in turn, is computed 
by multiplying the observed spectrum at the near detector (located at 300 m from the neutrino source) 
by the so-called 'far-near ratio', i?, defined as the ratio between the predicted flux at the far and near 
detectors. This factor corrects for the fact that at the near detector, the neutrino source is not point-like, 
but sensitive to effects such as the finite size of the decay tunnel, etc., whereas at the Super-Kamiokande 
site the neutrino source can be considered as point-like. According to the neutrino oscillation parameters 
measured in atmospheric neutrino experiments 0] the distortion of the spectrum measured with the far 
detector is predicted to be maximal in the energy range between 0.5 and 1 GeV. The determination of 
R is the leading energy-dependent systematic error in the K2K analysis [3 El • 

The HARP experiment has a large acceptance in the momentum and angular range relevant for K2K 
neutrino flux. It covers 80% of the total neutrino flux in the near detector and in the relevant region for 
neutrino oscillations. Thus, it can provide an independent, and more precise, measurement of the pion 
yield needed as input to the calculation of the K2K far-near ratio than that currently available. 

The neutrino beam of the K2K experiment originates from the decay of light hadrons, produced by 
exposing an aluminium target to a proton beam of momentum 12.9 GeV/c. In this paper, the mea- 
surement of the double-differential cross-section, tPa^^ /dpdfl of positive pion production for protons of 
12.9 GeV/c momentum impinging on a thin Al target of 5% nuclear interaction length (Ai) is presented, 
i.e. reproducing closely the conditions of the K2K beam- line for the production of secondaries. 

The HARP apparatus p ^ is a large-acceptance spectrometer consisting of a forward and large-angle 
detection system. The forward spectrometer covers polar angles up to 250 mrad which is well matched 
to the angular range of interest for the K2K beam line. 

The results reported here are based on data taken in 2002 in the T9 beam of the CERN PS. About 
3.4 million incoming protons were selected. After cuts, 209 929 secondary tracks reconstructed in the 
forward spectrometer were used in this analysis. The results are given in the region relevant for K2K, 
that is the momentum range from 0.75 GeV/c to 6.5 GeV/c and within an angular range from 30 mrad 
to 210 mrad. The absolute normalization was performed using 280 542 'minimum-bias' triggers. 

This paper is organized as follows. The experimental apparatus is outlined in Section |21 Section |31 
describes tracking with the forward spectrometer. Section0]discusses the calculation of the reconstruction 
efficiency. Section |S1 summarizes the particle identification (PID) capabilities of the spectrometer and 
describes the PID algorithm. Sections [HI and [7| give details of the cross-section calculation. Results 
are discussed in Section |S1 A comparison with previous data is presented in Section ^ An illustrative 
calculation of the K2K far-near ratio is shown in Sectional A summary is given in Section [TTl 



2 Experimental apparatus 

The HARP detector, shown in Fig. consists of forward and large-angle detection systems. The con- 
vention used for the coordinate system is also given in the Figure. In the large-angle region a TPC 
positioned in a solenoidal magnet is used for tracking. The forward spectrometer is built around a dipole 
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Figure 1: Schematic layout of the HARP spectrometer. The convention for the coordinate system is 
shown in the lower-right corner. The three most downstream (unlabelled) drift chambers are only partly 
equipped with electronics and not used for tracking. 

magnet with an integral field of J BydL=0.66 Tm for momentum analysis, with large planar drift cham- 
bers (NDC) for particle tracking, and three detectors used for particle identification: a time-of-flight 
wall (TOFW), a threshold Cherenkov detector (CHE), and an electromagnetic calorimeter (ECAL). The 
target itself is located inside the TPC. Beam instrumentation, including three timing detectors (BTOF) 
and two threshold Cherenkov detectors (BCA and BCB), provides identification of the incoming particle 
and the determination of the interaction time at the target. The impact point of the beam particle on the 
target and its direction are measured by a set of multi-wire proportional chambers (MWPCs). Several 
trigger detectors are available to select events with an interaction and to define the normalization. 

Data were taken with several beam momenta and target configurations. In addition to the data taken with 
the thin aluminium target of 5% Ai at an incident proton momentum of 12.9 GeV/c, runs were also taken 
with an empty target holder. These data allow a subtraction to be made of the interactions occurring 
in the material on the path of the incident beam. Other relevant configurations for the measurement 
described here are the data taken with and without target with other beam momenta (1.5, 3.0, 5.0, 8.0, 
8.9 and 15 GeV/c) with electrons, pions and protons. These settings have been used to determine the 
response of the spectrometer to these particles in terms of efficiency, momentum resolution and particle 
identification capability. Data with thick Al targets, such as a replica of the K2K target, have also been 
taken, but are not yet used in the present analysis. The momentum definition of the T9 beam is known 
with a precision of the order of 1% 0. 

A detailed description of the HARP experiment is given in Ref. [S]. In this analysis we utilize primarily 
the detector components of the forward spectrometer and the beam instrumentation. Below, the elements 
which are important for this analysis will be mentioned. 



2.1 Beam and trigger detectors 

A schematic picture of the equipment in the beam line is shown in Fig. |3 It is instrumented with the 
following systems: 

• A set of four multi-wire proportional chambers (MWPCs) measures the position and direction of 
the incoming beam particles, with an accuracy of wl mm in position and ^0.2 mrad in angle per 
projection. 
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Figure 2: Schematic view of the trigger and beam equipment. The description is given in the text. The 
beam enters from the right. The MWPCs are numbered: 1, 4, 2, 3 from right to left. On the left, the 
position of the target inside the inner field cage of the TPC is shown. 



• A beam time-of-flight system (BTOF) measures time difference over a distance of 21.4 m. It 
is made of two identical scintillation hodoscopes, TOFA and TOFB (originally built for the NA52 
experiment 0), which, together with a small target-defining trigger counter (TDS, also used for the 
trigger and described below), provide particle identification at low energies. This allows separation 
of pions, kaons and protons up to 5 GeV/c and provides the initial time at the interaction vertex 
(to). The resolution is shown in Fig.|21 The to, i.e. the time at which the incident beam particle is 
predicted to cross the mid-planc of the target {z = 0), is calculated after particle identification. The 
weighted average of the individual measurements of to from the three timing detectors is calculated, 
taking into account the velocity, /3, of the particle using the known beam momentum and the particle 
mass deduced after identification. The timing resolution of the combined BTOF system is about 
70 ps. 

• A system of two N2-filled Cherenkov detectors (BCA and BCB) is used to tag electrons at low 
energies and to tag pions at higher energies. The electron and pion tagging efficiency is found to 
be close to 100%. At momenta larger than 12 GeV/c it is also possible to tag kaons as can be 
seen in Fig. 0] which shows the pulse height spectrum of BCA and BCB for a 12.9 GeV/c beam. 
This spectrum displays raw channel counts without pedestal subtraction (the pedestal is around 
channel 110). The kaon and pion peaks can be clearly distinguished from the pedestal peak at 
low pulse-height which is due to heavier particles below Cherenkov threshold such as protons. The 
electrons are part of the pion peak. 



The target is positioned inside the inner field cage of the TPC. It has a cylindrical shape with a nominal 
diameter of 30 mm. The aluminium (99.999% pure) target used for the measurement described here has 
a nominal thickness of 5% Aj. Precise measurements of the thickness have been performed at different 
locations on its surface and show a maximum variation between 19.73 mm and 19.85 mm. 

A set of trigger detectors completes the beam instrumentation: a thin scintillator slab covering the full 
aperture of the last quadrupole magnet in the beam line to start the trigger logic decision (BS); a small 
scintillator disk, TDS mentioned above, positioned upstream of the target to ensure that only particles 
hitting the target cause a trigger; and 'halo' counters (scintillators with a hole to let the beam particles 
pass) to veto particles too far away from the beam axis. 

The TDS is designed to have a very high efficiency (measured to be 99.9%). It is located as near as 
possible to the entrance of the TPC and has a 20 mm diameter, smaller than the target which has a 
30 mm diameter. Its time resolution 130 ps) is sufficiently good to be used as an additional detector 
for the BTOF system. 

A double plane of scintillation counters (FTP), positioned upstream of the dipolc magnet, is used to 
select events with an interaction in the target and outgoing charged particles in the forward region. The 
plane covers the full aperture of the dipole magnet, with the exception of a central hole with a diameter 
of 60 mm to let the beam particles pass. The efficiency of the FTP was measured using events which had 
been taken simultaneously using triggers which did not require the FTP and amounts to >99.8%. 
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Figure 3: The timing resolution of the beam TOF detectors. The left hand panel shows the time difference 
measured between TOFA and TOFB, the right hand panel the time difference between TDS and TOFA. 



2.2 Drift chambers 

The main tracking device of the HARP forward spectrometer is a set of large drift chambers (NDC) 
placed upstream and downstream of the dipole magnet. These chambers were originally built for the 
NOMAD experiment [5], where they served both as a target for neutrino interactions and as a tracker 
for the produced charged particles. 

The spectrometer contains five NDC modules, each of which is made of four chambers. The chambers 
consist of three wire planes, with one plane (x) of wires oriented vertically; the other two (u and v) 
are rotated with respect to the vertical by ±5 degrees. They have been described elsewhere |Sj and we 
refer to |S] for a detailed description of their performance under HARP conditions. This performance 
can be summarized in terms of two quantities, spatial resolution and hit efficiency per plane. After 
internal alignment of the individual wires, the spatial resolution of the chamber is about 340 /im. The hit 
efficiency is smaller in HARP than it was in NOMAD due to the use of a different, non-flammable but 
less efiicient gas mixture. The hit efficiency varies between 80% and 85% in the central NDC modules. 



2.3 PID detectors 

Particle identification is performed in the forward spectrometer through the combination of several detec- 
tors downstream of the dipole magnet (CHE, TOFW and ECAL). We refer to ^ for a detailed description 
of the three systems. 

A large scintillator wall (TOFW) covering the full acceptance of the downstream tracking system is used 
in conjunction with the timing information from the beam detectors to measure the time-of-fiight of the 
secondary particles for momenta up to 5 GeV/c. The TOFW measures the time-of-flight of particles 
emanating from the target, and this, together with the charged track trajectory length, I, determines the 
velocity, (3, of the particle. 

The single scintillator counters are BC408 bars from Bicron, 2.5 cm thick and 21 cm wide. The counters 
are grouped into three mechanical structures (palisades). In the left and right palisades, scintillators 
are 250 cm long and are mounted vertically, while in the central palisade scintillators are 180 cm long 
and are mounted horizontally. The counters overlap partially by 2.5 cm to ensure full coverage. The 
scintillator slabs are viewed by one photo-multiplier tube (Philips XP2020) at each side. The TOFW and 
its performance is described in detail in Ref. [TO]. With an intrinsic resolution of the individual counters 
of 160 ps, a time-of-flight resolution better than 180 ps is achieved using this detector in combination 
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Beam Cerenkov A Beam Cerenkov B 

Figure 4: Pulse height spectra expressed in ADC counts from beam Cherenkov counters BCA (left) and 
BCB (right) at 12.9 GeV/c. The pedestal (marked 'protons' representing particles with a velocity below 
threshold) is around channel 110 in both cases. 



with the BTOF system. 

The threshold Cherenkov detector (CHE) separates pions from protons for momenta above the pion 
threshold (2.6 GeV/c) and identifies electrons below the pion threshold. The radiator gas (perfluorobutane 
C4F10) is chosen for its high refractive index, which allows the detector to be operated at atmospheric 
pressure. The particles traverse about 2 m of the radiating medium and generate photons that are 
deflected by about 135° upward or downward by two large cylindrical mirrors 6 m long with a radius of 
curvature of 2.4 m. Thirty-eight EMI 9356-KA photo-multipliers were used for their very low noise and 
high gain characteristics. In order to increase their useful light-collection area to a diameter of 340 mm, 
the photo-multipliers were matched to aluminized Winston cones. 

Finally, the electromagnetic calorimeter (ECAL) provides electron rejection. It is segmented longitudi- 
nally into two planes. The two calorimeter planes were assembled from existing calorimeter modules of 
the CHORUS experiment JI]. These planes consist of 62 and 80 modules, covering a total active width 
of 4.96 m and 6.4 m, respectively. Each module is composed of scintillating fibres (1 mm diameter) 
embedded in extruded lead sheets with a volume ratio 1/4. The ratio of the energy deposition in the two 
planes is different for electrons compared to hadrons. In addition, the comparison of the momentum of 
the particle measured by the curvature of its trajectory and the energy deposition in the calorimeter pro- 
vides another way to identify electrons. The ECAL complements the electron rejection of the Cherenkov 
above the pion Cherenkov threshold. 



3 Tracking with the forward spectrometer 
3.1 Tracking algorithm 

The track reconstruction algorithm starts by building two-dimensional (2D) segments per NDC module. 
Those are later combined to create 3D track segments (also per module). The requirements are the 
following: 

• Plane (2D) segment: At least three hits out of four in the same projection (u, x or v) compatible 
with being aligned. The drift sign associated to each hit is decided during the plane segment 
reconstruction phase. 
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• Track (3D) segment: Two or three plane segments of different projections, whose intersection defines 
a 3D straight fine. In the case where only two plane segments are found, an additional hit in the 
remaining projection is required. This hit must intersect the 3D straight line defined by the other 
two projections. 

Consequently, to form a track segment at least seven hits (from a total of 12 measurement planes) are 
needed within the same NDC module. Once track segments are formed in the individual modules they 
are combined (downstream of the dipole magnet) to obtain longer track segments. Finally, downstream 
tracks are connected with either the interaction vertex or a 3D track segment in NDCl (the NDC module 
upstream the dipole magnet, see Fig.^l to measure the momentum. All these tasks are performed by a 
sophisticated fitting, extrapolation and matching package called RecPack ^21; which is based on the well 
known Kalman Filter technique |13| . 

The interaction vertex in this analysis is well defined. The transverse coordinates {x, y) are obtained by 
extrapolating the trajectory of the incoming beam particle, measured with the MWPCs (with an error of 
the order of 1 mm) , and the z coordinate can be taken as that of the nominal plane of the target (which 
is 19.80 mm thick). 

Consequently, the momentum of a track can be determined by imposing the constraint that it emanates 
from the vertex, that is, by connecting a 3D segment downstream the dipole magnet with a 3D point 
upstream the magnet. Tracks of this type are called 'VERTEX2 tracks', and the estimator of the mo- 
mentum obtained by connecting a 3D segment with the vertex 3D point is denoted '^2'- Specifically, this 
is done by extrapolating the downstream 3D segment to the nominal plane of the target, and imposing 
that the distance between the transverse coordinates thus obtained {xs,ys) and the {x,y) coordinates 
defined above is less than 10 cm (in practice one builds a which also takes into account the measure- 
ment errors). Tracks which extrapolate to distances larger than 10 cm are not considered (in fact, the 
inefficiency of the p2 algorithm, a few percent comes almost exclusively from this source) . 

Alternatively, one can measure the momentum connecting a 3D segment downstream of the dipole with 
a 3D segment in the NDCl module. These are called 'VERTEX4 tracks', and the estimator of the 
momentum is denoted '^4'. The way a downstream segment is connected with a NDCl 3D segment is 
described with some detail below. In essence, one requires a reasonable collinearity in the non-bending 
plane and obtains the momentum from the curvature in the bending plane. 

The availability of two independent momentum estimators allows the tracking efficiency to be measured 
from the data themselves. This is possible, since, a) the reconstruction methods providing the estimators 
P2 and p4 are independent b) P2 and have a Gaussian distribution around the true momentum p 
(The distribution is expected to be Gaussian in the variable 1/p rather than p. With the relatively good 
resolution the difference is negligible.) This makes it possible to use one of the estimators (^2) to measure 
the yields while the other (^4) is used to measure tracking efficiency. The estimator p2 is preferred to 
measure yields since it does not involve the use of the NDCl module, where tracking efficiency is lower 
than in the downstream modules (see the discussion on tracking efficiency in Section!^. 

Figure [S] shows the inclusive p2 distribution measured for events with an FTP trigger. The hole in the 
FTP largely suppresses the peak of beam protons at 12.9 GeV/c. The remaining peak corresponds to 
events with an FTP trigger caused by elastically scattered protons, protons with multiple scattering in the 
tail of the angular distribution, and by protons accompanied by soft particles produced upstream of the 
FTP. The linear correlation between p2 and p4, shown in Fig. (left panel) for simulated tracks, illustrates 
the fact that both are estimators of the same quantity, while the correlation between p4 and p (Fig. |^ 
right panel) shows that both are unbiased estimators of p. The small non-linearities and disagreements 
between p2 and p^ and between p^ and p have negligible contribution to the total systematic uncertainty. 

3.2 Momentum and angular resolution 

Following the previous discussion, the momentum measurement for VERTEX2 tracks is performed by 
extrapolating tracks built downstream of the magnet to the vertex plane. The algorithm performs a 
loop over allowed momenta. For each value of p one computes the extrapolated position {xt,yt) at the 
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Figure 5: Inclusive p2 momentum distribution of reconstructed tracks (see text for definition). The peak 
from non-interacting beam particles is visible at 12.9 GeV/c. 

target reference Zo coordinate and the matching with the event (xq^uq) coordinates. The momentum 
is then calculated by minimizing this x^- For VERTEX4 the algorithm is similar, but in this case, a 3D 
segment downstream of the magnet is matched to a 3D segment in NDCl module. In both cases, an 
upper cut on the minimum matching decides whether the matching is accepted or not. This reduces 
the background from tertiary particles (not coming from the primary vertex) in the case of matching 
with the vertex (VERTEX2 tracks). For VERTEX4, this cut reduces the background from particles 
interacting in the region between NDCl and NDC2. 

The momentum resolution as a function of the momentum is shown in Fig. [7| (left panel) , for the case of 
P2- The resolution can be measured using beam particles of several momenta. Also shown (open circles) 
is the corresponding resolution found using the Monte Carlo. 

The momentum resolution does not improve below 3 GeV/c due to details of the momentum reconstruc- 
tion algorithm and also because the particles traverse the material at larger angles so that the multiple 
scattering term is no longer a constant. This feature is well reproduced by the simulation. Figure (right 
panel), shows the angular resolution. Both the momentum and angular resolutions are small compared 
with the size of the bins used in this analysis (500 MeV/c momentum bins, up to 4 GeV/c, 1000 MeV/c 
from 4 to 5 GeV/c, 1500 MeV/c from 5 to 6.5 GeV/c, and 30 mrad angular bins). In the region of 
interest, the agreement between data and Monte Carlo is good for the momentum resolution, while for 
the angular resolution the difference is less than 1 mrad, negligible compared to the bin size. Thus effects 
due to the finite resolution are small, and it is safe to apply a Monte Carlo based correction. 

The charge misidentification rate has been estimated by computing the fraction of protons that are 
reconstructed with negative charge. This is done by measuring the fraction of negative particles with 
momenta above the pion CHE threshold that give no signal in CHE. The upper limit of 0.5% for the 
charge misidentification probability is found to be consistent with the known CHE inefficiency 



3.3 Definition of kinematical variables 

The final cross section, being rotationally invariant around the beam axis, can be expressed in polar 
coordinates {p, 9), where p is the true total momentum of the particle and 9 is the true angle with respect 
to the beam axis (approximately equivalent to the z axis). However, given the rectangular geometry of 
the dipole and of the drift chambers, some of the corrections needed to compute the cross-section are 
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Figure 6: Left panel: The correlation between p2 and P4, showing that both are estimators of the same 
quantity; right panel: the correlation between and p shows that p^ is an unbiased estimator of the 
momentum p within the momentum resolution and binning (from 500 MeV/c up to 1.5 GeV/c) used in 
the analysis. 
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Figure 7: Left panel: momentum resolution (p2) obtained from fits to data (points with error bars) 
taken using several well-defined discrete beam momenta and no target. Also shown (open circles) is the 
corresponding resolution found using the Monte Carlo. Right panel: angular resolution obtained from fits 
to data (points with error bars) taken using several well-defined discrete beam momenta and no target. 
The open circles show again the corresponding resolution found using the Monte Carlo. In the region 
of interest, the agreement between data and Monte Carlo is good for momentum resolution with values 
much smaller than the binning used in the analysis (from 500 MeV/c for p < 4 GeV/c up to 1.5 GeV/c at 
p = 6.5 GeV/c). Similarly the difference between measured and predicted angular resolution is negligible 
compared to the 30 mrad binning adopted in the analysis (see text). 
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most naturally expressed in terms of {p, 9^, Oy), where Ox = a,Tcta,ii(j)x/pz) and 9y — arctan(py/pz). Thus 
the conversion from rectangular to polar coordinates is carried out at a later stage of the analysis. 



4 Track reconstruction efficiency 

The track reconstruction efficiency, e*''^'^'^^{p,9x,dy), is defined as the fraction of tracked particles (with 
position and momentum measured) iV'"'^'^'^ with respect to the total number of particles iVP^'''^ reaching 
the fiducial volume of the HARP spectrometer as a function of the true momentum, p, and angles, 9x,0y'- 

' ^P'^^'^^>~ NP^^^^{p,0x,0y) ■ 

The track reconstruction efhciency can be computed using the redundancy of the drift chambers taking 
advantage of the multiple techniques used for the track reconstruction. The rest of this section details 
the steps leading to this calculation. The efficiency was calculated for positively charged particles only. 



4.1 The use of the p4 estimator to measure tracking efficiencies 

The calculation of the cross-section requires the knowledge of tracking efficiency and acceptance in terms 
of the true kinematical variables of the particle. Strictly speaking, this is only possible if one uses the 
Monte Carlo to compute these quantities. This would make the calculation sensitive to the details of the 
Monte Carlo simulation of the spectrometer. 

The existence of two independent estimators of the momentum allows the tracking efficiency to be mea- 
sured in terms of p4, taking advantage of the fact that it is Gaussian distributed around p and therefore 
can be used to approximate the latter. 

Therefore, a sample of p4 tracks is selected with well measured momentum imposing the additional 
constraint that the tracks emanate from the primary vertex. This is achieved by requiring that the 
distance of the track extrapolation to the MWPC vertex is smaller than 10 mm. By construction, the 
vertex of VERTEX2 tracks clusters at a small radius around the nominal vertex origin (defined by the 
MWPC resolution) which is fully covered by these VERTEX4 tracks. 



4.2 Module efficiency 

Using the selected sample of VERTEX4 tracks, one can measure the tracking efficiency and acceptance 
of individual NDC modules in terms of VERTEX4 kinematical quantities. 

The measurement of p2 requires a downstream segment which is then connected to the event vertex. In 
turn, a downstream segment can be made of a segment in NDC2, or a segment in any of the modules 
downstream NDC2 (that is NDC3, NDC4 and NDC5, see Fig. |SJ). For the purpose of the analysis one 
can treat conceptually those three NDC modules as a single module which we call back-plane. Thus, 
a downstream segment is defined as a NDC2 segment, a back-plane segment or a long segment which 
combines both NDC2 and back-plane. In all cases the measurement of P2 requires that the true particle 
has crossed NDC2. By definition the control sample of VERTEX4 tracks verifies that the true particle 
crossed NDC2 (since the measurement of P4 requires a downstream segment connected to NDCl), but 
not necessarily that a segment was reconstructed in NDC2 (since a good VERTEX4 track can be built 
with a back-plane segment and a segment in NDCl). In practice, the required condition is that at least 
six hits are found inside the road defined by the extrapolation of the downstream track to NDC2. It was 
verified that this condition has a negligible effect on the efficiency determination. 

The NDC2 efficiency £2 is defined as the number of segments reconstructed in NDC2 (in terms of VER- 
TEX4 kinematical quantities, p, 6x,0y) divided by the number of tracks in the VERTEX4 control sample. 
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Figure 8: Schematic layout of the downstream modules of the drift chambers (top view). 
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Figure 9: Segment efficiency of NDC2 module (see text for definition), as a function of p (left panel), 9x 
(centre panel) and 9y (right panel). The efficiency is computed as a three-dimensional function of the 
above variables. The plots show the individual projections. Points with error bars correspond to data, 
the dashed line to Monte Carlo. The agreement between Monte Carlo and data calculation is good except 
in the region of large 9^ and small p (see text). 



This is equivalent to finding the number of segments reconstructed in NDC2 divided by the number of 
particles reaching NDC2. Thus £2 measures the 'true' tracking efficiency of the NDC2 module, unfolded 
from other effects such as acceptance, absorption or decay. If a particle decays or is absorbed before 
reaching NDC2 it will not be included in the control sample and therefore it will not be included in the 
calculation of £2- 

Figure |5| shows £2 as a function of p, 9^, 9y (estimated from the VERTEX4 control sample). As expected 
the distribution is flat in terms of all three variables. The dots represent the calculation from the data 
themselves, and show that the NDC2 tracking efficiency is essentially 100%. The dashed line represents 
the Monte Carlo calculation, which agrees with the data calculation except in the region of low momentum 
and large, positive 9x- The inefficiencies in these regions are correlated and are due to edge effects which 
are not perfectly described in the Monte Carlo. This region is not used in the analysis. 

The back-plane efficiency, £{,, is defined as the number of segments reconstructed in the back-plane divided 
by the number of tracks in the VERTEX4 control sample. This definition folds tracking efficiency 
with acceptance and other effects such as absorption or decay (for example, one could have a well 
reconstructed VERTEX4 track with a NDC2 segment and a NDCl segment, decaying or undergoing 
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Figure 10: Segment efficiency of the back-plane, defined by modules NDC3, NDC4 and NDC5 shown 
as a function of p (left panel) , Ox (center panel) and 9y (right panel) . The efficiency is computed as a 
three-dimensional function of the above variables. The plots show the individual projections. Points with 
error bars correspond to data, the dashed line to Monte Carlo. 



a nuclear interaction in NDC2). The Monte Carlo tends to overestimate the efficiency by less than 5% 
on average (see Fig. Il()|l . At large positive 9x the Monte Carlo predicts a rise in the efficiency which is 
not seen in the data. This is a region where large angle tracks are further deflected by the magnet and 
traverse the drift chambers at large angles. As will be shown below, owing to the redundancy of the 
chambers, the overall downstream efficiency is well reproduced by the simulation. 

The knowledge of the efficiency of NDCl is not needed, since the yields are computed in terms of p2 
which does not use it. However, one could choose to measure the yields in terms of pi and use p2 as an 
estimator of true momentum to measure tracking efficiency. In that case, NDCl would play the same 
role that NDC2 plays in the current approach. 

Indeed, it is illustrative to compute £i for VERTEX4 tracks in terms of VERTEX2 kinematical quantities, 
P2, O'x: 6y)- The efficiency ei is defined as the number of segments reconstructed in NDCl divided by the 
number of tracks in a VERTEX2 control sample. This is equivalent to requiring the number of segments 
reconstructed in NDCl over the number of particles reaching NDCl. Thus ei measures the 'true' tracking 
efficiency of the NDCl module, unfolded from other effects such as acceptance, absorption or decay. 

Figure ITTl shows £i as a function of p,9x,0y (estimated from the VERTEX2 control sample). The 
distribution is relatively flat in 6y, while it has a marked dependence on 9x, with a minimum at 9x — 0. 
The distribution of ei as a function of 9x shows the inefficiency of NDCl associated with a saturation of 
the chambers due to the primary beam intensity. NDCl is the only module affected by this saturation 
effect, which is negligible downstream of the dipole magnet (as proved by inspection of £2)- The Monte 
Carlo tends to overestimate the efficiency by 15%. The dip in the efficiency at 9x close to the origin is 
induced by the saturation effect of the beam. This feature is simulated in the Monte Carlo by artificialy 
lowering the efficiency of the drift regions most traversed by undeflected beam particles. Although one 
can measure the NDCl tracking efficiency with good precision using the data, the marked dependence 
on 9x (which translates also in a dependence on p) suggests as a better strategy to measure the yields in 
terms of VERTEX2, which does not use NDCl. 



4.3 Downstream tracking efficiency 

The downstream tracking efficiency, ^ is defined as the number of tracks reconstructed downstream 
the dipole magnet (those include NDC2 single segments, back-plane single segments and NDC2-back- 
plane combined segments) divided by the number of particles reaching NDC2 (which defines the fiducial 
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Figure 11: Segment efRciency of NDCl (see text for definition) as a function of p (left panel), 6^ (center 
panel) and 9y (right panel). The efRciency is computed as a three-dimensional function of the above 
variables. The plots show the individual projections. Points with error bars correspond to data, the 
dashed line to Monte Carlo. 

volume). Since £2 and £{, are uncorrelated, this quantity is easily computed from the individual segment 
efficiency described above, as: 



_down 



£2 



£2 ■ £b 



(2) 



Figure 1121 shows the total downstream segment efficiency in terms of the kinematical quantities (from 
VERTEX4) p, 9x,9y. Due to the very high tracking efficiency of the individual modules (the apparent drop 
of the efficiency of the back-plane visible in Fig.^|is largely due to decay and absorption) the downstream 
tracking efficiency is almost 100%, flat in all variables, within the angular acceptance considered. The 
overall downstream tracking efficiency is well reproduced by the Monte Carlo. 

4.4 Upstream tracking efRciency 

For the contribution of the vertex matching to the overall track reconstruction efficiency, one needs to 
compute the fraction of times that a good downstream track segment was correctly matched to the vertex 
point resulting in an acceptable momentum. This is also done using a VERTEX4 control sample, with 
the additional constraint that the particle emanates from the target volume. The latter condition is 
ensured by requiring rv4 < 30 mm, x^atch — where rv4 is the distance of the track extrapolation to 
the MWPC vertex in the target reference plane, and Xmatch the goodness of the matching of the track 
with the beam particle track extrapolated to the nominal target position. The cut in xf^^tch guarantees 
good matching with the vertex. These tracks enter the denominator of the efficiency calculation. The 
numerator is made of the accepted VERTEX2 tracks within this sample: 



^(3p2;3p4,rv4 < 30 mm,xLtch < 10) 



^vertex4(gp^^^y^ < 30 mm, xLtch < 10) 



(3) 



where the condition 3p4(2) guarantees that the track momentum was estimated by the two reconstruction 
algorithms, respectively. 

The efficiency for reconstructing a VERTEX2 track for a particle coming from the target is shown in 
Fig. ^] Figure [TSl (upper-left), shows that the efficiency as a function of is flat and close to 95%, 
up to 150 mrad and drops above this value. This drop is due a the momentum-dependent acceptance 
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Figure 12: Downstream tracking efficiency as a function of kinematic variables, p, 0^, and 6y. at produc- 
tion for positively charged particles emanating from the vertex. Upper- left panel: As a function of 9x- 
Upper-right panel: As a function of 6y. Lower- left panel: As a function of p. Lower-right panel: As a 
function of p averaged over the 9x and 9y regions used in the present analysis only. The efficiency is flat 
in all variables and close to 100%. The solid histograms correspond to data, the dashed line to Monte 
Carlo. 
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Figure 13: Upstream tracking efRciency as a function of kinematic variables, p, 9^, and 9y, at production 
for positively charged particles emanating from the vertex. Upper-left panel: As a function of O^- Upper- 
right panel: As a function of By. Lower- left panel: As a function of p. Lower-right panel: As a function 
of p averaged over the and dy regions used in the present analysis only. The efficiency is close to 100% 
for negative and momenta less than 4 GeV/c, and drops for high values of 6^ (due to the dipole's 
acceptance) and for high momenta (due to a weakness in the reconstruction algorithm, see text). Points 
with error bars correspond to data, the dashed line to Monte Carlo. The agreement is excellent, except 
in the momentum bin 6.5 GeV/c < p < 8.0 GeV/c, where there is a 25% difference. This bin is not used 
in the analysis (see text). 



limitation imposed by the dipole magnet, as clearly demonstrated by Fig. ^] (upper-right) , which shows 
a flat distribution in the non-bending plane 0y and Fig. 1131 flower- left V which shows the efRciency as a 
function of p, integrated for all 9x- The drop of the efficiency for large values of 6x is fully correlated 
with the drop at low p. This can be seen in Fig. 1131 flower- right), where the reconstruction efRciency as a 
function of p for negative 9x (particles fully contained in the dipole acceptance) is shown. The efRciency 
is flat and close to 100% up to 4 GeV/c, and drops for the last two bins. The efRciency as a function of 
9x {9y) for momenta less than 4 GeV/c is close to 100%, indicating that the loss of efRciency is due to the 
efRciency drop at p > 4 GeV/c. The efRciency of the momentum reconstruction algorithm is lower near 
the edge of the acceptance because it requires a number of trajectories around the best Rt to be inside 
the aperture. This requirement reflects in a drop of overall efflciency for 9x > 150 mrad. 

The drop in reconstruction efficiency at high momenta is due to the lack of optimization, at high mo- 
mentum, of the reconstruction algorithm. Its effect on the calculation of the cross-section is very small 
(in practice it translates into a slightly larger error in the higher momentum bins, which are dominated 
by statistical errors). Except for the highest-momentum bin, the data are well described by the Monte 
Carlo. 
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4.5 Total reconstruction efficiency 



The total tracking efBciency, s^'^'^'^^^ can be expressed as the product of two factors. One factor represents 
the downstream (of the dipole magnet) tracking efficiency and the other represents the efficiency for 
matching a downstream segment to a vertex (for tracks originating inside the target volume): 

Ty^down ;\rvcrtcx2 

^track ^down ^VGrtcx2 /^\ 

^parts ^down ' ^ ^ 

where iV™''*'''^^ is the number of VERTEX2 tracks, which corresponds to N^'^'^'^^ in Eq.^ In the momentum 
range of interest for this analysis, a good time-of- flight measurement is essential for particle identification. 
Therefore, one also requires a good TOFW hit matched to the track. In addition to selecting the 
scintillator slab hit by the particle, this matching consists of a cut in the matching of the track with 
the TOFW hit coordinate that is measured by the time difference of the signals at the two sides of the 
scintillators. 

The TOFW hit matching efficiency was computed using the same track sample as in the previous section: 

^ToF ^ jV^°^(3ToF;3p4,rv4 < 30 mm, xLtch < 10) .5^ 
^vcrtex4(3p^^ rv4 < 30 mm, xiatch < 10) ' ^ ' 

where the condition 3ToF guarantees that a TOFW hit was associated with the track, and the total 
reconstruction efficiency is found from: 

The total reconstruction efficiency is shown in Fig. El The inclusion of the TOF wall enhances the 
momentum-dependent acceptance cut (due to the fact that the TOF wall has a smaller geometrical 
acceptance than the NDC back-plane). The total reconstruction efficiency as a function of 9^ has a 
slight slope and drops for positive 9x above 100 mrad. The total reconstruction efficiency as a function 
of p for negative 6^ is flat (and about 90%) for momenta below 4 GeV/c. The drop to 90% is due to 
the inefficiency in matching tracks to a TOFW hit. The agreement between data and Monte Carlo for 
the total reconstruction efficiency is excellent, except in the last momentum bin, where there is a 25% 
difference, reflecting the same effect already observed in the upstream efficiency. 



5 Particle identification 



A set of efficient FID algorithms to select pions and reject other particles is required for the current 
analysis. A Monte Carlo prediction of the differential yields of the various particle types shows that the 
pion production cross-section is small above 6.5 GeV/c, which is set as the upper limit of this analysis. 
The electron distribution peaks at low energy, while the proton background increases with momentum. 
The kaon yield is expected to be only a small fraction of the pion yield. In the momentum and angular 
range covered by the present measurements the proton yield is of a similar order of magnitude as the 
pion yield. 

The PID strategy is based on the expectation of the yields of different particle types predicted by the 
Monte Carlo, and also on the momentum regions covered by the available PID detectors. The time-of- 
flight measurement with the combination of BTOF and TOFW systems (referred to as the TOFW mea- 
surement in what follows) allows pion-kaon and pion-proton separation to be performed up to 3 GeV/c 
and beyond 5 GeV/c respectively. The Cherenkov is used for hadron-electron separation below 2.5 GeV/c 
and pion-proton/kaon separation above 2.5 GeV/c in conjunction with the TOFW. The ECAL is used 
only to separate hadrons from electrons below 2.5 GeV/c to study the Cherenkov performance. 

As mentioned above, the electron (positron) background is concentrated at low momentum (p < 2.5 GeV/c) 
It can be suppressed to negligible level with an upper limit on the CHE signal, given the fact that electrons 
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Figure 14: Total reconstruction efficiency as a function of kinematic variables, p, 9x, and 9y, at production 
for positively charged particles emanating from the vertex. Upper- left panel: as a function of 9x- Upper- 
right panel: as a function of 9y. Lower- left panel: as a function of p. Lower- right panel: as a function of 
p averaged over the 9x and 9y regions used in the present analysis only. The efficiency is close to 90% 
for negative 0^ and momenta less than 4 GcV/c, and drops for high values of 9^ (due to the TOFW 
acceptance) and for high momenta (due to a weakness in the reconstruction algorithm, see text). Points 
with error bars correspond to data, the dashed line to Monte Carlo. The agreement is excellent, except 
in the bin with highest momentum, where the difference is 6%. 
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are the only particles giving signal in the Cherenkov below the pion Cherenkov light emission threshold, 
which is equal to 2.6 GeV/c for the gas mixture used in HARP. In practice, any particle that has a 
momentum below 2.5 GeV/c and a signal in the CHE exceeding 15 photo-electrons is called an electron. 
In the following we will refer to this cut as the e-veto cut. The remaining electron background after the 
e-veto cut is negligible as studied in Rcf. |14| . 

Having applied the e-veto cut to reject electrons and keeping in mind that there is a small fraction of 
kaons, one builds PID estimators for protons and pions by combining the information from TOFW and 
CHE using likelihood techniques fSec. \b.'2]i . Then, a cut on these PID estimators is applied to select 
pions or protons. The selected samples (raw pion and proton samples) will contain a small fraction of 
kaons, which can be estimated from the data, as described in Ref. TT. This background is subtracted 
from the dominant yields of pions and protons. 

The quantities that enter the cross-section calculation are the raw pion and proton yields and the PID 
efficiencies and purities (PID corrections) obtained by the application of the e-veto cut and cuts in the 
PID estimators. The PID corrections include the e-veto efficiencies, the kaon subtraction corrections and 
the pion-proton efficiency matrix (M*"*, described in Sec. 16.4(1 . A precise knowledge of these quantities 
(as a function of momentum and angle) requires the understanding of the responses of different PID 
detectors to the particle types considered. This is studied in detail using the data, taking advantage of 
the redimdancy between the PID detectors. 

These steps are explained briefly in the following sections. More details are given in Ref. |14j . 
5.1 Response of the PID detectors 

The TOFW-CHE ProbabiUty Density Function (PDF), P(/3, iVpho|«,P, 6*), describes the probability that 
a particle of type i (pion or proton) with momentum p and polar angle 9 results in simultaneous measure- 
ments f3 in TOFW and iVphc in CHE. The latter comes directly from the calibrated CHE signal, while (3 is 
the particle velocity, computed as /3 = ^tof/(itof 'c), where Itot is the track length measured from the nomi- 
nal vertex position to the TOFW hit position, itof is the measured time-of- flight and c is the speed of light. 
Assuming that the PID from both detectors are independent^ the TOFW-CHE PDFs can be factorized 
in independent TOFW and CHE PDFs, such that P(/3, 7Vphe|i,P, ^) = PiP\i,p,0) ■ P{Npi,^\i,p,9). 

5.1.1 TOFW response 

The use of the particle velocity, /3, to characterize the TOFW response has several advantages. Its 
distribution is nearly Gaussian-the agreement between data and Monte Carlo is sufficient for the current 
analysis-and it discriminates very effectively between pions and protons up to momenta around 5 GeV/c 
(at this energy the separation between the average values of the proton and pion Gaussians is around 
2.2f7). These points are illustrated in Fig. 1151 

To build the TOFW PDFs one should know the (3 distribution of pions and protons as a function of the 
particle momentum and angle. In order to maximize the efficiency of the selection algorithm and to avoid 
any possible bias in the PID corrections related with a data-MC disagreement, those distributions have 
been measured from the data. This requires the ability to select pure and unbiased samples of pions and 
protons from the data. Samples of pions with negligible contamination from other species can be obtained 
selecting particles of negative charge passing the e-veto cut (Fig.^J left panel). At low momentum, the 
proton parameters can be obtained by simply fitting to a Gaussian the proton part of the (5 distribution, 
since this is well separated from pions (Fig. ^1 central panel). At large momentum (> 2.5 GeV/c), the 
proton and pion distributions overlap significantly. In this case, pions are rejected to the 1% level by a 
cut on the CHE signal (Sec. I5.1.2|l . The resulting sample contains a majority of protons and residual 
contaminations from pions and kaons. In this momentum range the proton parameters can be obtained 

'^The use of the reconstructed momentum instead of the true momentum in the probabihty density function 
P(/9, AfphoK) Pi 6) introduces a small correlation between TOFW and CHE. However, the good momentum resolution of 
the forward spectrometer {dp/p < 10%) makes this correlation very small. At the level of precision required by this analysis 
this effect can be neglected. 
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Figure 15: The distribution of the particle velocity, /3, for pions and protons of 3 and 5 GeV/c, for 
both data (corresponding to a sample of beam particles, selected as pions or protons by the beam 
instrumentation) and its corresponding Monte Carlo simulation. At 3 GcV/c the separation between the 
two populations is ~ 5(T, and the separation is still « 2.2a at 5 GeV/c. The proton and pion peaks in 
the Monte Carlo have been separately normalized to the area of the corresponding peaks in the data. 



by fitting the inclusive /3 distribution to a triple Gaussian with fixed pion and kaon shapes^ , as shown in 
the right panel of Fig. ^| A full description of this technique is given in Ref. ^31 ■ 

Figure [T7I shows the mean value and the standard deviation of (3 for pions and protons. The left panel 
of Fig. El shows the comparison between the values obtained with the method explained above and the 
ones obtained using the MC information about the true particle type. In both cases, the same sample 
of Monte Carlo data has been used. The observed good agreement confirms that the proposed technique 
does not bias the pion and proton parameters. The result of the application of this method to the data 
is shown in the right panel of Fig. El This plot also shows the comparison between data and Monte 
Carlo when using the same method for the selection of pure particle samples. We observe a global shift of 
about 0.003 in /? which can be attributed either to a TOFW misalignment of 3 cm along the z direction 
or to a time offset of 0.1 ns between the TOFW and the BTOF system. In order to account for this 
difference the analysis uses the PDFs based on the f3 distribution measured in the data directly when 
treating 'data-events' and the PDFs based on the Monte Carlo distribution when treating 'Monte Carlo 
events'. The TOFW PDFs are parametrized as the sum of a dominant Gaussian function and a term 
accounting for non-Gaussian outliers, normalized to the observed effect. 



5.1.2 Cherenkov response 



The Cherenkov detector is used digitally in this analysis: a signal is accepted if the number of photoelec- 
trons is larger than 2. To obtain the CHE PDFs for a given momentum and angular bin, the fraction 
of true pions (protons) with negative (positive) signal in the CHE is measured. The number of true 
particles of a given type is obtained with a technique similar to the one used for the TOFW detector, in 
this case by applying a strict cut to the TOFW measurement (see Ref. [T^ for the details). Fig. 1181 fleft 
panel) shows the CHE inefficiency for particles of negative charge (essentially pions) as a function of the 
reconstructed momentum and angle. The asymptotic inefficiency for pions, i.e. above a momentum of 
3.5 GeV/c, is estimated to be (1.0 ± 0.5)%. 

In the momentum range studied no signal is expected in the CHE for protons. However, in a fraction of 
events, the reconstruction algorithm wrongly associates the CHE hit from a pion or an electron to the 

■^The pion shape is obtained from the data, while the kaon shape is calculated using the MC information. 
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Figure 16: Inclusive /? distribution for pions, kaons and protons passing the e-veto cut. Left panel: 
negative particles (essentially pions) with reconstructed momentum between 3.25 and 4 GeV/c. Central 
panel: positive particles with reconstructed momentum between 1.75 and 2.25 GeV/c. Right panel: 
positive particles (mostly protons) with reconstructed momentum between 3.25 and 4 GeV/c. In the 
right panel, the kaon and pion Gaussians are also shown (shaded areas). 
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Figure 17: Mean value of P (points) and its standard deviation (error bars and rectangles) as a function 
of momentum for protons and pions. On the left, MC events only, obtained using only reconstructed 
quantities to select pure samples (points with error bars) and using the MC information about the true 
particle type (rectangles). On the right, points with error bars correspond to similarly selected events 
from real data, and the rectangles correspond to equivalent MC events. The centres of the rectangles are 
indicated with points. 
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Figure 18: Left panel: CHE pion inefRciency as a function of the momentum for different angular regions. 
Right panel: CHE proton efficiency as a function of the angle for different momentum regions. 

proton and consequently, a fraction of protons has a non-negligible amount of associated photoelectrons. 
This is a potential source of background (as well as of pion inefficiency), particularly important at high 
momentum, where the TOFW is not applicable. The efficiency of the CHE for protons has been measured 
as a function of momentum and angle. The results are shown in Fig. 1181 fright panel). This non-zero 
efficiency is fully taken into account in this analysis as explained in Ref. 14 . The CHE PDFs are given 
in Fig. IIHI 



5.2 The pion-proton PID estimators 

The assignment of a particle type z = tt, p to a reconstructed track is based on a cut Pi > Pent in a PID 
estimator Pi, the so-called combined PID probability, which is built by using Bayes' theorem: 

P^PC\^N 0)^ P{/3,N^i.,\i,p,9)-P{t,p,e) 

n n»IP,^Vphc,P,(^J P(/3,iVphek,P,e)-P(^,P,e)+P(/3,7Vphc|p,P,e)-P(p,P,e) ' ^ ' 

where P{(3, A^phcKiP, are the TOFW-CHE PDFs described above, and P{i,p, 9) is the prior, describing 
the a priori probability that a particle passing the event and track selection criteria fSecs. l^Tl and l^?^ is 
of type i and has momentum p and polar angle 6. 

Several simplifications have been made to the general formula JTJ. They result in a slightly less efficient 
particle type selection, which implies a larger error on the PID corrections. However, the additional error 
is negligible when added in quadrature to the dominant non-PID errors. As a first approximation, the 
PDFs that enter the PID estimator were averaged over all angles. As second approximation, equally 
probable priors were used, so that they cancel. In this case only the information from the current track is 
used to build the PID estimator. The PID estimator built with no priors does not have a full probabilistic 
meaning and cannot be used directly to estimate the particle yields. Instead, the raw pion and proton 
yields must be corrected by the efficiencies and purities obtained by the application of the cut Pi > Pcut, 
as will be described in Sec. [731 Finally, the TOFW-CHE PDFs can be factorized in independent TOFW 
and CHE PDFs, as explained before. The final PID estimator is then represented by the formula: 

p_pnflA7 ^- P(/3|z,p)-P(jVpheN,p) 

^ P(/3k,p)-P(^phekp)+P(/3|p,p)-P(iVphe|p,p) ■ ^ ' 
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6 Calculation of the cross-section 



The double-differential cross-section for the production of a particle of type a can be expressed in the 
laboratory system as: 



where is expressed in bins of true momentum {pi), angle {0j) and particle type (a), and the terms 

on the right-hand side of the equation are: 

• N^j, is the number of particles of observed type a' in bins of reconstructed momentum (p^/) and 
angle {0j'). These particles must satisfy the event, track and PID selection criteria, explained below. 
This is the so called 'raw yield'. 

• M~jl.^,j,^, is a correction matrix which corrects for finite efficiency and resolution of the detector. It 
unfolds the true variables ija from the reconstructed variables i'j'a' and corrects the observed num- 
ber of particles to take into account effects such as reconstruction efficiency, acceptance, absorption, 
pion decay, tertiary production, PID efficiency and PID misidentification rate. 

• f/^p^ is the inverse of the number of target nuclei per unit area {A is the atomic mass, Na is the 
Avogadro number, p and t are the target density and thickness). 

• -^pot is the number of incident protons on target. 

The summation over reconstructed indices i'j'a' is implied in the equation. It should be noted that the 
experimental procedure bins the result initially in terms of the angular variable 9, while the final result 
will be expressed in terms of the solid angle O. Since the background from misidentified protons in the 
pion sample is not negligible, the pion and proton raw yields (N^j, , for a' — tt, p) have to be measured 
simultaneously. 

For practical reasons, the background due to interactions of the primary proton outside the target (called 
'Empty target background') has been taken out of the correction matrix M~^. Instead, a subtraction 
term is introduced in Eq. |^ 



dV„ _ 1 A 



7V^,,(T)-7Vf?,,(E) , (10) 



where (T) refers to the data taken with the aluminium target and (E) refers to the data taken with no 
target (Empty target). 

The event, track and particle identification selection criteria will be described first, then the method used 
to obtain the cross-section and each of the corrections will be described in more detail. 



6.1 Event Selection 

In the 12.9 GeV/c beam protons are selected by vetoing particles which give a signal in any of the beam 
Cherenkov detectors. Only particles which give a good timing signal in all three beam timing detectors, 
leave a single track in the MWPCs, and are not seen in the halo detectors are accepted. A good timing 
measurement is defined as a set of three hits, one in each of the timing detectors, with their relative time 
difference consistent with a beam particle. The distribution of the position of beam particles extrapolated 
to the target is shown in Fig. ^] (left panel). The size of the target is indicated by a circle. Only particles 
extrapolated within a radius of 10 mm are accepted. By evaluating the number of tracks reconstructed 
in the spectrometer as a function of the extrapolated impact point of the MWPC track to the target, it 
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Figure 19: Left panel: reconstructed position {xo,yo) of beam particles at the reference zq plane of the 
target. The circle gives the position and size of the target. Right panel: observed number of tracks per 
incident proton as a function of the maximum accepted measured radius of incidence of the beam track. 
The ratio is normalized to unity at 6 mm. 



Data Set Al 5% 12.9 GeV/c 12.9 GeV/c Empty Target 

Protons on target 17,954,688 4,769,408 

Total events processed 4,710,609 771,330 

Events with accepted beam proton 3,404,372 547,838 

Prescaled triggers with accepted beam proton 280,542 74,522 

FTP triggers 2,087,732 225,639 

FTP trigger rate = ( FTP triggers / pot ) 0.116 0.047 

Total good tracks 209,929 11,704 

Table 1: Total number of events in the 12.9 GeV/c aluminium 5% Ai target and empty target data sets, 
and the number of protons on target as calculated from the prescaled trigger count. 



was determined that (1.5 ±0.5)% of the proton tracks selected according to these criteria miss the target, 
as shown in Fig. 1191 fright panel). A correction for this loss has been applied. The MWPC track was 
required to have a measured direction within 5 mrad of the nominal beam direction to further reduce 
halo particles. The purity of this proton sample is estimated to be better than 99.5%. 

Prior to those cuts, the beam particle was required to satisfy the trigger conditions described in Sec- 
tion Applying the above selection cuts to the 12.9 GeV/c aluminium 5% Ai target data set and 
the 12.9 GeV/c empty target data set results in the total statistics listed in Tabled The total number 
of protons on target {Npot in Eq. ^ listed in the table is exactly the number to be used in the overall 
normalization of the cross-section results, and is known to better than 1%. The total number of protons 
on target is counted using prescaled 'beam' triggers that were continuously recorded at the time of data 
taking. The trigger condition for the prescaled beam triggers only involved a simple coincidence of scin- 
tillators in the beam line with no requirement of an interaction in the target. Using subsamples of the 
triggers the prescale factor was checked to confirm it had its preset value 1/64. Because the selection 
criteria for beam protons used in event analysis and prescaled beam proton events are the same, the 
efficiencies for these cuts cancel, and the total normalization can be known without additional systematic 
uncertainty. 

Events to be used in the analysis must also contain one or more hits in the forward trigger plane (FTP). 
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6.2 Track Selection 



The recorded events have been processed according to the track selection criteria hsted below: 

• The VERTEX2 track momentum is measured (see Section IS^ . 

• A track segment in NDC2 or in the back-plane is used in track reconstruction. 

• Number of hits in a road around the track in NDCl > 4 (this is applied to reduce non-target 
interaction backgrounds). 

• The average for hits with respect to the track in NDCl < 30, 

• Number of hits in the road around the track in NDC2 > 6 (this is applied to reduce background of 
tracks not coming from the target). 

• The track has a matched TOFW hit. 

The result of applying these cuts to the entire 12.9 GeV/c aluminium 5% Ai and empty target data sets 
is listed in Table 

In addition, geometrical cuts are applied. As described in Sec.0]for positive 6^ the efficiency is momentum 
dependent. This region is avoided in the analysis by defining the fiducial volume as —210 < 6*2; < mrad 
(thus, only particles in the negative half of the bending plane of the dipole are accepted) and —80 < 
9y < 80 mrad. The restricted acceptance in 9y is imposed to avoid edge effects of the dipole, possible 
fringe effects in the magnetic field, etc. Since the behaviour of the spectrometer is calibrated with beam 
particles (at 9y = 0) the analysis restricts 9y to a rather small region around the horizontal mid-plane 
of the spectrometer. In order to avoid a correction for the acceptance of the FTP-trigger and to avoid 
background from beam protons, cross-sections are given for 9 > SO mrad. 

6.3 PID selection 

Particle identification criteria (described in Sec. [SJ are applied to the tracks passing the event and track 
selection criteria. First the e-veto cut is applied to reject electrons and then a cut in the PID estimator is 
applied to distinguish between pions and protons. Figure 1201 shows the combined pion probability (PID 
estimator for pions) for positive particles passing the e-veto cut. A large population of particles in the 
low probability region is attributed to a contribution from protons. The small peaks at 0.5 correspond 
to particles which leave no useful information in either the TOFW or the Cherenkov. The peak near 0.9 
in the right panel corresponds to particles which leave no useful information in the TOFW (presumably 
being non-Gaussian outliers) but which give a positive signal in the CHE. 

It is found that the optimal cut to select pions with high efficiency and purity is Pj^ > 0.6. The cut is set 
at a value of probability where the track population is low, and thus the result is not sensitive to small 
changes in the exact value. Protons are selected by the condition P^r < 0.4 (equivalent to Pp > 0.6). 

6.4 The Atlantic and UFO analyses 

Two complementary analyses have been performed with the aim of checking internal consistency, and 
checking for possible biases in the respective procedures. The first, called Atlantic'^, simplifies the prob- 
lem of unfolding by decomposing the correction matrix of Eq. ((SJ into distinct independent contributions, 
which are computed mostly using the data themselves. The second analysis, called UFO (from UnFOld- 
ing), performs a simultaneous unfolding of p, 9 and PID, with a correction matrix computed mainly 
using the Monte Carlo. 

^'Atlantic' for Analysis of Tracks at Low ANgle with Tof Id and Cherenkov id 
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Figure 20: Combined pion probability (logarithmic scale) for positive particles passing the e-veto cut. 
The left panel shows positive particles below CHE threshold. The right panel shows positive particles 
above CHE threshold. The meaning of the different spikes is explained in the text. 



The UFO procedure uses an iterative Bayesian technique, described in Ref. Ts', in order to unfold the 
measured distribution. The central assumption of the method is that the probability density fmiction 
in the physical parameters ('physical distribution') can be approximated by a histogram with bins of 
sufficiently small width. A population in the physical distribution of events in a given cell ija generates a 
distribution in the measured variables, Mijai'j'a', where the indices ija indicate the binning in the physi- 
cal angular, momentum and PID variables, respectively, and i'j'a' the binning in the measured variables. 
Thus the observed distribution in the measurements can be represented by a linear superposition of such 
populations. The task of the unfolding procedure consists then in finding the number of events in the 
physical bins for which the predicted superposition in the measurement space gives the best description 
of the data. 

In order to predict the population of the migration matrix element Mijai'j'a': the resolution, efRciency 
and acceptance of the detector are obtained from the Monte Carlo. This is a reasonable approach, since 
the Monte Carlo simulation describes most of these quantities correctly (see Section 0J). Where some 
deviations from the control samples measured from the data are found, the data are used to introduce 
(small) corrections to the Monte Carlo. 

Although some corrections are common to both approaches, large differences between the results of these 
two analyses would indicate inconsistencies in the simplifications adopted by Atlantic for unfolding, the 
hypothesis of correct Monte Carlo description of the detector on which UFO is based, or both. As it 
turns out, the analyses are consistent within the overall systematic error, reinforcing our confidence in the 
correctness of the results presented here. For clarity, in the rest of this paper only the Atlantic analysis 
will be discussed. 



7 The Atlantic analysis 

As discussed in section 13.21 both the momentum and angular resolution are small compared with the 
binning of the cross-section. Migration effects are, therefore, small. In particular, angular migration can 
be neglected. In addition, kinematic migration is almost decoupled from pion-proton PID migration. As 
explained in Sec. |31 electron and kaon ID has been decoupled from the dominant pion-proton ID so that 
electron and kaon correction factors are diagonal in the PID variables. With the above considerations 
the correction matrix can be written as: 

^a.'/.' = (M^ao')-' • ■ (M^r)-' ■ iMl,)-' , (11) 
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where again reconstructed indices are indicated with a prime. The corrections are apphed in the order 
from right to left as they appear in the equation. The symbols in Eq. have the following meaning: 



^ija' the collection of factors applying the corrections that are diagonal in the PID indices: 
reconstruction efficiency, acceptance, physical loss of particles (absorption, decay), background 
from tertiary interactions, e-veto efficiency and kaon subtraction; 

{Mf^,)~^ is the simplified unfolding matrix correcting for the momentum smearing which only 
depends of the indices i and i' representing the true and reconstructed momentum bins, respectively; 

{Mjj,)~^ is the identity matrix, representing the assumption that the smearing effect in the angular 
measurement is negligible; and 

i-^ij-aa')^^ is the matrix which corrects for pion-proton PID inefficiency and migration, which is 
diagonal in i, i' and j, j', but built of two-by-two sub-matrices, each different and non-diagonal in 
the PID variables a, a'. 

The diagonal efficiency correction 

-1 _ rccon acc ^ absorption tcrtiarios , K o ('-|9N 

'^ija> ~ "^ij ^ij ^ija' "^ija' 'lija' 'lija' y'^'^) 

is composed of the following factors: 

^rocon correction for the overall reconstruction efficiency; 
wfj'^ the correction for the acceptance; 

^ifa°'''''°" the correction for the loss of particles due to absorption and decay; 

the correction for the background of tertiary particles generated by the secondaries pro- 
duced in the target; 

vfja' t^*^ factor correcting for the kaon background; and 
Vija' is i-^^ factor correcting for the effects of the electron veto. 

The first two corrections are the same for pions and protons while the latter four also depend on the 
particle type. It is worth noting that the efficiency correction is expressed in terms of the true momentum 
and angle, and in terms of the reconstructed particle type (a'). This is because these corrections are 
applied before PID unfolding, as explained below. 

As advanced in Sec. 13.31 some of the above corrections are computed as a function of p, 9^ and 6y, while 
some others are directly expressed in the final variables {p,0). In the first case, the transformation to 
polar coordinates (p, 9) is done integrating over all Ox and 9y resulting in a given 9 bin. In particular, 
the four first corrections of Eq. (|12|l . denoted by w, are computed as a function of (p, 9.j;, 9y). 

Each of the above corrections will be described in the sections below. 



7.1 Reconstruction and acceptance corrections 

The correction for the total reconstruction efficiency, requiring a momentum measured and a matched 
TOFW hit (computed in SecEj), is introduced as a weight wj^^™" = [e"™"(p, 6':r, 

It is necessary to correct for the restricted definition of fiducial volume. Inside the 9y acceptance (that 
is, below the vertical cutoff at 9y = ±80 mrad) the correction is a simple factor of 2 due to the fact that 
tracks with 9x > Q are not used. For values of 9 above the 9y cutoff, the correction is: 

1 /tan(6l,';"')\ 

s-^(p,..,.,) = -.arcsin(^-L|^) , (13) 
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Figure 21: Sketch of the forward detector showing the two kinds of acceptance corrections, a) for 9 < 6\ 
where e'^'^'^ = 0.5, and b) for 9 > 9y^*' where e'*'^'^ is given by Eq. H13|l. 



describing the part of the circle which is inside the acceptance. Figure |2] shows a sketch depicting the 
two forms of geometrical acceptance and the origins of the correction factors listed above. The acceptance 
correction is then applied as a weight defined in Eq. (|12|l . w^j'^ — [e^'^'^{p, 9^, 9y)]~'^. The above corrections 
are independent of the particle type. 



7.2 Corrections for absorption, decay, and secondary interactions 

An additional correction to consider is the absorption and decay of secondary pions (protons) in the 
materials of the detector components upstream of the magnet which prevent them from reaching the 
downstream region of the detector. These missing particles are not considered by the reconstruction 
efficiency described above. In some cases, e.g. pion decay, tertiaries are reconstructed as part of the 
original track. A correction is applied to take these cases into account. The overall absorption and decay 
rate is determined using the Monte Carlo, as a three-dimensional function e'^^°'^^^^°'^{p,9x,9y), whose 
projections are shown in Fig. [53 The overall effect is between 10% and 30% depending on p and 9x. 
This effect is verified to be correct within 10% of its magnitude using beam particles. The 'absorption' 
correction is then applied as a weight introduced in Eq. = [1 - e'^^°'P^'°'^{p, 9x,9y)]'^ for 

both particle types separately. 

A correction of opposite sign to the one above stems from positively charged particles which are not 
produced in the primary interaction between the incident proton and the target nucleus. These ter- 
tiary particles ('tertiaries') can come directly from the target (nuclear re-interactions, which is a small 
effect, since the target is only 5% Ai in thickness) or from the region outside the target area. Another 
different background is due to particles produced by the interactions of the primary proton with ma- 
terial outside the target; this is corrected by taking data with empty target settings, and is described 
in Section [731 The correction for tertiaries is also computed with Monte Carlo, as a three-dimensional 
function etei-tiarics^p^ Q^^^ shown in Fig. [23 The overall effect is between 2-3% (for pions) and 7% (for 
protons). In addition one needs to correct for pion decay resulting in muons which tend to be coUinear 
with the original pion. The 'tertiary' correction is then applied as a weight introduced in Eq. H12|l . 
Wija^'^"'^^ = 1 — e|^°''*"*'''°^(p, 9x,9y). The correction is model-dependent, and has been assigned a system- 
atic uncertainty of 100%. 

The above corrections are computed separately for true pions and protons. However, for practical reasons, 
they are applied before FID unfolding assuming that they correspond to reconstructed pions and protons 
(hence the index a' in Eq. I12|l . The bias introduced by this approximation is negligible since pion-proton 
mixing is very small ( < 5%), as demonstrated in Fig. 1261 and these corrections are either small (< 7% for 
tertiaries) or similar for both particle types (the absorption is similar while the decay of pions introduces 
a relatively small correction). Thus, the maximum bias would be of the order of 0.4%. 
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Figure 22: Absorption and decay rate as a function of kinematic variables, p, 9^, and 9y, at production 
for positively charged particles emanating from the vertex. Solid lines are pions, dashed lines are protons. 
Left panel: As a function of p. Central panel: as a function of 6x- Right panel: As a function of 6y. The 
acceptance effects for positive 9x due to particles hitting the dipole walls are clearly visible. This effect 
is momentum dependent, as seen in the left panel. 
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Figure 23: Tertiary particle rate as a function of kinematic variables, p, 6^, and 6y, at production for 

positively charged particles emanating from the vertex. Circles are pions, squares are protons, triangles 
are muons falsely identified as pions. Left panel: As a function of p. Central panel: as a function of ^x- 
Right panel: As a function of 9y. 
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7.3 Empty target subtraction 



Several additional background sources need to be corrected for. The 'empty target' background is defined 
as the particles accepted by the selection criteria which are generated by interactions of the primary pro- 
tons outside the target. The effect of this background is measured experimentally to a good approximation 
by taking data without placing the target in its holder. 

The corrections described in the previous section are applied to both 5% Ai target and empty target data 
sets. The empty target yield N°tj, (E) undergoes similar corrections to the yields measured with the target 
in place. The corrected empty target yields are then subtracted bin-by-bin from the corrected yields 
measured with target to remove this background. The relative normalization of the data with target and 
the empty target data is calculated using the number of protons on target accepted in the prescaled beam 
trigger. The overall subtraction is approximately 20% as shown in Fig. 1241 The approximation used in 
this approach is to assume that the target itself does not influence the primary proton beam. To first 
order, this assumption introduces an error of 5% on the subtraction, given by the interaction length of 
the target. 



^ 0.002 

o 






> 0.0018 






O 0.0016 




30 - 60 mrad 


0.0014 




— Target 


O 0.0012 




Empty 


a 




0.001 




Subtracted 


*^ 0.0008 






0.0006 






0.0004 






0.0002 


1 

r 1 .... 1 .... 1 .... 1 ... . \-v;-\-v\-v; 



0.002 
0.0018 
0.0016 
0.0014 
0.0012 

0.001 
0.0008 
0.0006 
0.0004 
0.0002 



60 - 90 mrad 

- - - Target 

Empty 

Subtracted 




momentum (GeV/c) 



momentum (GeV/c) 



0.002 
0.0018 
0.0016 
0.0014 
0.0012 

0.001 
0.0008 
0.0006 
0.0004 
0.0002 









90-120 mrad 




— Target 




Empty 




Subtracted 








' 1 



^ 0.002 

o 

> 0.0018 
O 

O 0.0016 

^ 0.0014 

O 0.0012 
Q. 

^ 0.001 
^ 0.0008 
0.0006 
0.0004 
0.0002 









120-150 mrad 




— Target 




Empty 




— Subtracted 




1 



momentum (GeV/c) 



momentum (GeV/c) 



o 

> 0.0018 
O 0.0016 




150-180 mrad 


0.0014 




— Target 


9 0.0012 

a 




Empty 


~~- 0.001 




; Subtracted 


^ 0.0008 






0.0006 






0.0004 






0.0002 






■ 1 ■ ■ ■ ■ 1 ■ ■ ■ ''"\":^•!•\••.••^ ...I-.-.- r : n-.-.- 



Q t ■■rvrt ...I ..l..i_J 

1 2 3 4 5 6 

momentum (GeV/c) 



^ 0.002^ 

o 

> 0.0018 '- 




momentum (GeV/c) 



Figure 24: Positive pion yields, defined as pions per p.o.t (protons on target) as a function of momentum. 
Empty target yields are subtracted from target yields to remove backgrounds. The subtraction ranges 
from negligible to approximately 20%. 
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7.4 Corrections for electron-veto and kaon background 



As discussed, the electron and kaon hypotheses are not considered by the PID selection algorithm. Kaons 
produced in the target and identified as either pious or protons have to be subtracted. Electrons are 
rejected by applying the CHE veto described above. This veto introduces a loss of efficiency for pious 
and protons. Multiplicative corrections are applied to the raw yields of pious and protons to compensate 
for these two effects as shown in Fig. 1251 

The electron veto, described in Sectional has some effect on pion and proton efficiency (around 7-10% 
of the pious and protons also give a signal in CHE below pion threshold, due to spatially associated 
electrons, e.g. the emission of hard (5-rays) which has been measured using both the Monte Carlo and the 
data (pure electron and hadron control samples were selected for that purpose using ECAL The 
weights -qlj^, (Eq. (fTT|l ) are the inverse of the efficiencies for pious and protons to survive the electron 
veto requirement, respectively. 

Kaons are subtracted from the raw yields of pious and protons. The kaon yield has being estimated from 
the data by fitting the inclusive T0FW-/3 distribution to the sum of three Gaussians (corresponding 
to protons, kaons and pious). An example is shown in Fig. 1161 The kaon-to-pion and kaon-to-proton 
migration rates have also been determined from the data using the same technique as for pion-proton 
migration. The subtraction can be expressed as a correction factor rj^^,, introduced in Eq. (|ll|l . The 
correction factors for pious and protons are shown in the bottom panels of Fig. 1251 The correction for 
pious is only relevant (~ 3%) in the region of transition between TOFW and CHE (^ 3GeV/c), while 
the correction for protons is of the order of 10% in the entire phase space. 

A detailed description of the methods used to determine the electron and kaon correction factors can be 
found in |14j . 



7.5 Unfolding the momentum dependence and the pion— proton yields 

A yield of tracks with reconstructed momentum p2 falling in a bin i' can be expressed as a superposition 
of tracks with true momentum in a bin i [pi). The coefficients of this expansion are the elements of the 
momentum migration matrix of Eq. 111|) . Mf^,, that is n{p2)i' = Mf^, ■ n{p)i. To perform the momentum 
unfolding each observed track with measured quantities {p2, 0x, Oy) populates several bins in a histogram 
of true variables {p, 9) with weights M^^, . This is mathematically equivalent to the matrix inversion of 
Eq. fTT)l . 

The raw yield of identified pious and protons once corrected for all terms that are diagonal in the PID 
variables, nf- = e~^^, ■ (Mf^, )-i • (M^^^, • [7V"^, (T) - {E)\ , is related to the true pion and proton yield, 
n^-, by the PID migration matrix (also called PID efficiency matrix), M^j.^^,, introduced in Sec. 16.41 In 
each bin of p and 9 the matrix M^'^ is defined by: 

where the elements of the matrix, in rows, are the fractions of observed pious that are true pious (Mtt^), 
observed pious that are true protons (M^rp), observed protons that are true pious (Mp^r), and observed 
protons that are true protons (Afpp). In this experiment this matrix can be computed using the redun- 
dancy in the data, as described in il4i. Figure shows the elements of Af' (and the corresponding 
errors), as a function of momentum, for different angular intervals. 

Then, the true yields can be computed by solving the system of linear equations given by Eq. H14() . The 
covariance matrix of the true yield vector, n", is computed by error propagation, taking into account the 
covariance matrices of AP"^ and of the observed yield vector, n°' : 
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Figure 25: On the top, pion (left) and proton (right) c-vcto correction factors averaged over angles. On 
the bottom the kaon-pion (left) and kaon-proton (right) correction factors in real data for ^ = 90 mrad. 



33 



1.05 
1 

0.95 
0.9 

0.85 
0.8 







fEt+pt 



13 = 45 mrad 

I I ^3= 1 50 mrad 

I I I I I I I I I I I I I 



P (GeV/c) 




EEE 

i.J^ 



_L 



2 4 6 8 

P (GeV/c) 



+ 



P (GeV/c) 



P (GeV/c) 



Figure 26: The elements of the PID efheiency matrix in the data, as a function of momentum, for two 
different angles. The probability cut is placed at 0.6. 
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Figure 27: Measurement of the double-differential 7r+ production cross-section in the laboratory system 
d?a/{dpdfl) for incoming protons of 12.9 GeV/c on an aluminium target as a function of pion momentum 
p, in bins of pion polar angle 9. The data points are the measurements, the histogram represents the 
Sanford-Wang parametrization fitted to the data. 



where all indices run over the pion and proton hypotheses. 



8 Results 



Figure [23 and Table [3 show the measurement of the double-differential cross-section for positive pion 
production in the laboratory system as a function of the momentum and the polar angle. Only diagonal 
errors are shown in the plots and table (a full discussion of the error evaluation is given below). Also 
shown in Figure|23is a fit to a Sanford-Wang parametrization, which will also be discussed in this section. 
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Table 2: HARP results for the double-differential tt"*" production cross-section in the laboratory system, 
(fa"" /{dpdn). Each row refers to a different (Pmin ^ P < Pmax, ^min ^ ^ < ^max) bin, where p and 9 
arc the pion momentum and polar angle, respectively. The central values quoted are the ones obtained 
via the Atlantic analysis discussed in the text. The square-root of the diagonal elements of the covariance 
matrix are also given. 
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8.1 Error estimates 



A detailed error analysis has been performed to evaluate the accuracy of the pion cross-section measure- 
ment. The main errors entering in this measurement are listed below. 

First, the statistical uncertainties associated with the track yields measured from the aluminium target 
setting and from the empty target setting (needed for subtraction, as explained above) have been included 
in the pion production cross-section uncertainty estimates. 

Second, several uncertainties associated with the corrections needed to convert the measured track yields 
to true track yields have been evaluated. The track reconstruction efficiency correction is based on the 
combination of thin target aluminium and beryllium data sets. The main error associated with this 
computation is given by the size of the statistical sample. The correction to the pion and proton yields 
due to absorption or decay is computed via a Monte Carlo simulation. An uncertainty of 10% for both 
proton and pion yields has been assumed for this correction, in addition to the uncertainty due to the 
finite size of the simulated data sample used to estimate this correction. Similarly, simulated data (and 
their associated uncertainties) were used to estimate the correction for the contamination in the sample 
due to tertiary particles that are not produced in the target, but rather by the decay of secondaries, or 
by the interaction of secondaries in the spectrometer material. An uncertainty of 100% has been assumed 
for this subtraction, for both proton and pion yields. Furthermore, an uncertainty has been assigned to 
the empty target subtraction, in order to account for the effect of the target itself which attenuates the 
proton beam. 

Third, uncertainties associated with the particle identification of tracks, and with the corrections needed 
to convert yields of tracks identified as pious to true pion yields, have been included. Among the several 
error sources associated with the pion-proton FID selection, the dominant one is due to the uncertainty 
in the (small) fraction of pions and protons with an associated anomalous TOFW /3 measurement, that 
is a /3 measurement which exhibits a non-Gaussian behaviour. Estimates of the uncertainty in the kaon 
contamination and in the correction for the electron veto have been obtained from an analysis of the data 
as explained in Section|Sl The robustness of the pion FID selection and its associated correction has been 
evaluated by performing the analysis with tighter and looser FID probability cuts with respect to their 
nominal values, while correcting for the FID efficiency and migration corresponding to the probability 
cuts. 

Fourth, we have included uncertainties associated with the momentum reconstruction performance of 
the spectrometer, and with the corrections needed to convert the measured momenta to 'true' momenta. 
Concerning the momentum, biases and resolution effects are taken into account using both real and 
simulated data. It was found that momentum biases do not exceed the 5% level from a study of beam 
particles at different momenta and from a comparison between the reconstructed momenta and the 
momenta inferred from (3 measurements with the TOFW and the threshold curves in the Cherenkov. 

Finally, an overall normalization uncertainty of 4% has been estimated. The dominant sources for this 
uncertainty are the targeting efSciency uncertainty, which is deduced from the measurement of transverse 
beam spot size on target, as well as the reconstruction and FID uncertainties that are fully correlated 
across different {jp, 9) pion bins, and which are not included in the above evaluation. On the other 
hand, the aluminium target thickness and density were carefully measured, and the effect on the overall 
cross-section normalization due to these uncertainties is negligible. 

8.2 Results of the Error Evaluation 

The impact of the error sources discussed in the previous section on the final cross-section measurement 
has been evaluated, either by analytic error propagation, or by Monte Carlo techniques. Correlation 
effects among different particle types, and among different (p, 9) bins, have also been taken into account. 

The cross-section uncertainty level is quantified by adopting two different conventions. The rationale 
is that both the errors on the 'point-to-point', double-differential cross-section, and the error on the 
cross-section integrated over the entire pion phase space measured, might be of interest. 
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Error Category 


Error Source 


^diff (%) 


r ('07' ^ 


Statistical 


Al target statistics 


1.6 


0.3 




bmpty target subtraction (stat.j 


1.3 


0.2 




Sub-total 


2.1 


0.4 


Track yield corrections 


Reconstruction efficiency 


0.8 


0.4 




Pion, proton absorption 


2.4 


2.6 




Tertiary subtraction 


3.2 


2.9 




bmpty target subtraction (syst.) 


1.2 


1.1 




Sub-total 


4.5 


4.1 


Particle identification 


PID Probability cut 


0.2 


0.2 




Kaon subtraction 


0.3 


0.1 




Electron veto 


2.1 


0.5 




Pion, proton ID correction 


2.5 


0.4 




Sub-total 


3.5 


0.7 


Momentum reconstruction 


Momentum scale 


3.0 


0.3 




Momentum resolution 


0.6 


0.6 




Sub-total 


3.2 


0.7 


Overall normalization 


Sub-total 


4.0 


4.0 


All 


Total 


8.2 


5.8 



Table 3: Summary of the uncertainties affecting the double-differential cross-section ((^cjiff) and integrated 
cross-section (^jnt) measurements. See text for details. 



First, the dimensionless quantity (5(j;g- is defined, expressing the typical error on the double-differential 
cross-section, as follows: 

where i labels a given pion {p, 9) bin, {A'^a^ / {dp-dVl))i is the central value for the double-differential cross- 
section measurement in that bin, and ((5[A^(t'^/ {dp ■ d^)])i is the error associated with this measurement. 



The individual and cumulative effect of the error sources discussed above on the 5^]g quantity are shown 
in Table O The typical error on the double-differential cross-section is about 8.2%. The dominant error 
contributions to 5^^^ arise from overall normalization (4%), subtraction of tertiary tracks (3.2%), and 
momentum scale (3.0%). More details on the relative double-differential cross-section uncertainties are 
shown in Fig. |2H1 for all measured {Pt9) bins. In Figure ESI and in Tab. |21 the individual cross-section 
uncertainties are grouped into five categories: statistical, track yield corrections, particle identification, 
momentum reconstruction, and overall normalization uncertainties. Uncertainties associated with the 
track yield corrections discussed above dominate the cross-section uncertainties in the low momentum 
region, while the dominant errors in the high momentum region are due to the momentum reconstruction 
and to the overall normalization. 

Second, we define the dimensionless quantity 5^^^, expressing the fractional error on the integrated pion 
cross-section, cr'^ (0.75 GeV/c <p< 6.5 GeV/c, 30 mrad < < 210 mrad), as follows: 

_ y^E..,(ApA17),C,,(ApA17), 

where {A'^a^)i is the double-differential cross-section in bin i, (A^cr'^/(ApAr2))i, multiplied by its cor- 
responding phase space element (ApAri)^. Here, Cij is the covariance matrix of the double-differential 
cross-section obtained by summing thirteen matrices from the error sources listed in Table O and whose 
square root of the diagonal elements, ^/Cu, corresponds to the error {d{A'^a^ / {ApAil)))i appearing in 
Eg. 1151 This covariance matrix is used to compare the two independent analyses of the same cross-section 
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Figure 28: Estimate of the fractional errors on the double-differential pion production cross-section 
measured as a function of pion momentum p and polar angle. The errors shown are in percent. The 
contributions from the error categories in Tab. |21 (thin histograms) as well as the total error (thick solid) 
are shown. 
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Parameter Value 

(4.4 ± 1.3) • 10^ 

C2 (8.5 ± 3.4) • 10-1 

C3 (5.1 ±1.3) 

C4 = C5 (1.78 ±0.75) 

C6 (4.43 ±0.31) 

C7 (1.35 ± 0.29) • 10-1 

C8 (3.57 ± 0.96) • 10^ 

Table 4: Sanford-Wang parameters and errors obtained by fitting the dataset. The errors refer to the 
68.27% confidence level for seven parameters (Ax^ — 8.18). 



measurement, and to obtain the best-fit values, errors, and correlations for the coefficients entering into 
the Sanford-Wang formula used to parametrize the HARP measurements. The correlation coefficients 
among distinct {p,9) bins in Cy vary between —0.19 and +0.95. 

The contributions to (5jjj^ from all the error sources considered, as well as the total error estimate on the 
integrated cross-section, are also given in TableOl As expected, (mostly) correlated errors such as the one 
from the normalization or tertiary subtraction remain (almost) as large as they were for the point-to-point 
error. On the other hand, the contribution of the momentum scale uncertainty is negligible here, since 
its effect tends to be anti-correlated among different phase space bins. In addition to the normalization 
and tertiary subtraction, other uncertainty sources which have some impact on the integrated cross- 
section include the pion absorption correction and the empty target subtraction systematic uncertainty. 
Overall, the total uncertainty on the pion production cross-section measured over the entire phase space 
(0.75 < p < 6.5 GeV/c, 30 < < 210 mrad) is estimated to be about 6%. 

In the following section, the cross-section results are also expressed in a parametrized form. 



8.3 Sanford-Wang parametrization 

The TT^ production data was fitted with a Sanford-Wang parametrization |17) , which has the functional 
form: 

^ -1-77^ -KP^ ^) = (^IP ) exp[-C3-J5 C66'(p - CTPbeam COS= S)] , (17) 

where X denotes any system of other particles in the final state, Pbeam is the proton beam momentum in 
GeV/c, p and 9 are the tt"*" momentum and angle in units of GeV/c and radians, respectively, cPa / (dpdfl) 
is expressed in units of mb/(GeV/c sr), dfl = 2i: d(cos6), and the parameters ci, . . . ,C8 are obtained 
from fits to 7r+ production data. 

The parameter ci is an overall normalization factor, the four parameters 02,03,04,05 can be interpreted 
as describing the momentum distribution of the secondary pions, and the three parameters cg,C7,C8 as 
describing the angular distribution for fixed secondary and proton beam momenta, p and Pbeam- 

This 

formula is purely empirical. In the minimization procedure, seven out of these eight parameters were 
allowed to vary. The parameter 05 was fixed to the conventional value C5 = 04, since the cross-section 
dependence on the proton beam momentum cannot be addressed by the present HARP data-set, which 
includes exclusively measurements taken at Pbeam — 12.9 GeV/o. In the minimization, the full error 
matrix was used. 

Concerning the Sanford-Wang parameters estimation, the best-fit values of the Sanford-Wang parameters 
are reported in Tabled together with their errors. The fit parameter errors are estimated by requiring 
Ax^ = x^~Xmin~^-l^' corresponding to the 68.27% confidence level region for seven variable parameters. 
Significant correlations among fit parameters are found, as shown by the correlation matrix given in 
Table El 
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Parameter 


Cl 


C2 


C3 


C4 = C5 


C6 


C7 


C8 


Cl 


1.000 














C2 


-0.056 


1.000 












C3 


-0.145 


-0.691 


1.000 










C4 = C5 


-0.322 


-0.890 


0.831 


1.000 








C6 


-0.347 


0.263 


-0.252 


-0.067 


1.000 






C7 


-0.740 


0.148 


-0.067 


0.077 


0.326 


1.000 




C8 


0.130 


-0.044 


0.205 


-0.040 


-0.650 


0.189 


1.000 



Table 5: Correlation coefficients among the Sanford-Wang parameters, obtained by fitting the data. 
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Figure 29: Projections of the differential light hadron production cross-section as a function of p inte- 
grated over the range 30 < < 210 mrad (left panel), and production cross-section as a function of 6 in 
the range 0.75 < p < 6.5 GeV/c (right panel). The points show the HARP measurements, the dotted 
curve the best-fit Sanford-Wang parametrization. 



The HARP cross-section measurement is compared to the best-fit Sanford-Wang parametrization of 
Table m in Figs. EH and EHl 

The goodness-of-fit of the Sanford-Wang parametrization hypothesis for the HARP results can be assessed 
by considering the best-fit value of x^jn=305 for 41 degrees of freedom, indicating a very poor fit 
quality. We note that the goodness-of-fit strongly depends on the correlations among the HARP cross- 
section uncertainties in different (p, 9) bins, and therefore cannot be inferred from Fig. 1271 alone. If these 
uncertainties were (incorrectly) treated as completely uncorrelated, the best-fit value would decrease 
from 305 to 57. A more comprehensive study of tt"*" production at various beam momenta and from 
various nuclear targets in HARP is planned and will follow in a subsequent publication, and should 
hopefully shed more light on the cause of the poor quality of the Sanford-Wang hypothesis reported here. 



9 Comparison with existing forward pion production data on 
aluminium 

Finally the HARP results are compared with existing 7r+ production data available in the literature 
directly from aluminium targets ^1 El QOl 1^ ■ The comparison is restricted to proton beam momenta 
between 10 and 15 GeV/c (close to the K2K beam momentum of 12.9 GeV/c), and for pion polar angles 
below 200 mrad (the range measured by HARP and of relevance to K2K). 

The comparison is based on the HARP Sanford-Wang parametrization rather than on the HARP data 
points themselves, in order to match pion momenta and angles measured in past Al experiments. Further- 
more, a correction to rescale the HARP Sanford-Wang parametrization at 12.9 GeV/c beam momentum 
to the 10-15 GeV/c beam momenta of the past Al datasets is applied |22| . 
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Figure 30: Comparison of the double-differential pion production cross-section measured in HARP, and 
the one measured in past experiments using an aluminium target and 10-15 GeV/c momentum beam 
protons. The points are the data from past experiments, and the shaded area reflect their normalization 
uncertainty. The solid line is the HARP Sanford-Wang parametrization rescaled to the beam momentum 
of past experiments, as discussed in the text. 



Given these model-dependent corrections, it was found that the HARP results are consistent with Ref. 
and Refs. 2t|, agree rather well with ,21, and are somewhat lower than, but still marginally consistent 
with, Ref. . Figure QUI shows the comparison between HARP and the above datasets. 



10 HARP results as input to the K2K far— to— near neutrino flux 
ratio prediction 

The main application of the measurement presented in this paper, the double-differential 7r+ production 
cross-section, (Pa (p -I- Al ^ 7r+ + X)/ (dpdil), is to predict the far-near ratio, R, for the muon neutrino 
disappearance search in the K2K experiment. 

As discussed in Sec.^ the determination of the far-near ratio is the leading energy-dependent systematic 
error in the K2K analysis [21 13 • To compute this quantity a Monte Carlo program simulating all relevant 
beam-line geometry and materials, and all relevant physics processes, is used. In this simulation, the 
neutrino flux prediction uncertainty is dominated by the uncertainties in the forward tt^ production 
arising from the interactions of the 12.9 GeV/c protons in the aluminium target material. Therefore, it 
is instructive to recompute a prediction for the K2K far-to-near flux ratio prediction based on the new 
experimental information presented in this paper. The HARP-based prediction has been obtained by 
substituting the original 7r+ production cross-section assumed in the K2K beam Monte Carlo hadronic 
model with the HARP Sanford-Wang parametrization discussed in Section [SJ while keeping unchanged 
all other ingredients of the K2K beam MC simulation, such as primary beam optics, pion re-interactions 
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in the aluminium target, pion focusing, pion decay, etc. More details on the default K2K beam MC 
assumptions can be found in Ref. [2]. 




J^tGeV) E,(GeV) 



Figure 31: Muon neutrino fluxes in the K2K experiment as a function of neutrino energy iJ^, as predicted 
by the default hadronic model assumptions in the K2K beam Monte Carlo simulation (dotted histograms), 
and by the HARP 7r+ production measurement (full circles with error bars). Left panel shows unit-area 
normalized flux predictions at the K2K near (top) and far (bottom) detector locations, $near and $fg^j., 
respectively, while right panel shows the far-to-near flux ratio ^fg^j./'l'near (open squares with error boxes 
show K2K model results). 

The result of this exercise is shown in Fig. |^ The left panel shows muon neutrino fluxes in the K2K 
experiment as a function of neutrino energy E^, as predicted by the default hadronic model assump- 
tions in the K2K beam Monte Carlo simulation (dotted histograms), and by the HARP 7r+ production 
measurement (filled circles). The plots on that panel show unit-area normalized flux predictions at the 
K2K near (top) and far (bottom) detector locations, <I>near and ^f^r, respectively. Right panel shows the 
far-to-near flux ratio $fj^j./$near obtained from K2K Monte Carlo (empty squares with error boxes) and 
from HARP measurament (filled dots with error bars). The fluxes predicted by HARP and the present 
K2K model are in good agreement within the errors. This is reflected also in a good agreement in R, 
in particular in the oscillation region (below 1.5 GeV). Finally, it is worth noting that the error on R 
associated with the HARP measurement (including statistical and systematic errors) is of the order of 
1%, since most errors on the cross-section cancel in the ratio. The current systematic error attached to R 
in the K2K analysis is of the order of 7%. Thus, although the result presented here does not yet represent 
a new measurement of R (which requires a full evaluation of other systematic errors independent of the 
HARP measurement but associated with the K2K beam line setup), it clearly shows the considerable 
improvement that can be achieved by K2K by using this new measurement. In addition, the data taken 
with the replica of the K2K target will be valuable to study the effect of reinteractions which needs to 
be taken into account in the beam simulation. 



11 Summary and conclusions 

In this paper we present a measurement of the double-differential production cross-section cPa^^ /dpdfl, 
for positively charged pions. The incident particles are protons of 12.9 GeV/c momentum hitting a 
thin aluminium target of 5% Aj. The measurement of this cross-section has a direct application to the 
calculation of the neutrino flux of the K2K experiment. The data were taken in 2002 in the T9 beam of the 
CERN PS. Out of 4.7 million triggers processed, 3.4 million incoming protons were selected. After cuts. 
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around 210 000 secondary tracks reconstructed in the forward spectrometer were used in this analysis. 
These high statistics results were corrected for measurement resolutions. These data were fitted with a 
Sanford-Wang parametrization. The results are given for positively charged pions within a momentum 
range from 0.75 GeV/c to 6.5 GeV/c, and within an angular range from 30 mrad to 210 mrad. The 
average statistical error is 1.6% per point. The absolute normalization was performed using prescaled 
beam triggers. The overall efficiency for track reconstruction and particle identification is known to better 
than 6%, while the average point-to-point error is 8.2%. 
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